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inflation ~ 60 e-foldings is predicted: this is an improvement over classical measures that generally 
are ambiguous as to whether inflation is certain to occur. 
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I. INTRODUCTION 

Quantum cosmology is regarded as a possible way of 
obtaining the initial conditions required to start the evo- 
lution of a classical cosmological model, for a general in- 
troduction, see Ref. m. These initial conditions corre- 
spond to a number of arbitrary constants that determine 
entirely, in the absence of chaos, the future evolution of 
the model. In the simple Friedman-Lemaitre-Robertson- 
Walker (FLRW) metric these constants determine the 
amount of matter present, and in the case of a scalar field 
source the initial balance of kinetic energy to potential 
energy of the scalar field. This determines how the initial 
expansion proceeds and whether the strong-energy con- 
dition is first violated to create an inflationary expansion. 
One could hope to overcome ambiguities in the predic- 
tion of whether inflation occurs that result from purely 
classical measures of probability 

The spatial curvature k is also another constant to be 
supplied. Although with strong-energy satisfying matter 
e.g. radiation, it will only dominate at large scale fac- 
tor to either force the universe to re-collapse or to keep 
eternally expanding. But when the strong-energy condi- 
tion is violated the curvature instead dominates at small 
scale factors. Curvature is therefore especially important 
if the universe is assumed to be created in an inflationary 
state. But, because closed k = 1 models have finite size 
they have generally been thought to be most relevant for 
quantum cosmology as they will have finite action, see 
e.g. Ref. |. 

The archetypal model that has been considered is the 
closed DeSitter space that has at small scale factors a 



Euclidean or "forbidden region" . Starting at zero scale 
factor one can envision tunneling into the classically al- 
lowed region. Boundary conditions need to be supplied in 
order to make predictions: the two most widely used be- 
ing the Hartle-Hawking (HH) jij and Tunneling ones || . 
The Tunneling boundary condition can be formulated in 
a number of ways, see [|6[ for a recent discussion, but we 
will generally refer to pj where the Tunneling condition is 
defined as "outgoing modes" only. Roughly speaking the 
HH condition makes use of the Euclidean nature of space 
time to smooth out singularities while the Tunneling one 
corresponds to allowing only outgoing modes, making it 
analogous to quantum a decay of an atom. 

To overcome the limitation of only having closed mod- 
els quantum creation of bubbles during any subsequent 
inflationary phase can enable locally open regions to form 
. The use of these so-called Hawking- Turok instantons 
are essentially making use of the fact that in DeSit- 
ter space all curvatures are equivalent. Different slicings 
of the 5-dimensional DeSitter hyperboloid correspond to 
different curvature k when considered as a 4-dimensional 
model, see e.g. Ref. ]io|| . But if the final requirement is 
simply an open or flat universe this would seem, at best, 
a rather convoluted procedure for their production. Be- 
cause open universes can also be compact [[fl) it seems 
possible to work directly with a universe of specific cur- 
vature, instead of by necessity starting with a closed uni- 
verse that can later by quantum tunneling create locally 
open regions. This will be the subject of this paper, for a 
general reference to topology in cosmology see For 
some results this requirement of compactness seems un- 
necessary as one can let the volume become arbitrarily 
large. Or alternately quantize the Friedmann equation 
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for the scale factor directly where the volume factor is a 
redundant multiplicative factor cf. [ Bj] . 

Allowing k = — 1 has a drastic effect on quantum cre- 
ation scenarios as the forbidden region that is assumed 
to be tunneled through is no longer present. The clas- 
sical singularity at zero scale factor is no longer isolated 
from the larger universe but instead classical evolution 
can start from arbitrarily small size. 

The fact that curvature plays such an important role 
in quantum creation scenarios is worrying since the stan- 
dard notions of curved space, and their associated met- 
rics, are likely to require extensive modifications as the 
quantum gravity epoch is approached. As a first approx- 
imation to including the effect of quantization of the cur- 
vature, we will also consider that the curvature is initially 
a quantum (q) variable that is allowed to take values 
within an ensemble. With this assumption we will argue 
that creation is in some sense "more likely" if this initial 
tunneling is not required which leads one to conclude that 
open universes are favoured. A somewhat related work 
p3[ has concluded that quantum tunneling is favoured if 
the strong-energy condition is only just being violated. 
This can be understood since this corresponds to making 
the barrier to be tunneled through very shallow and so 
easily overcome. We are taking this idea a step further 
by removing the barrier entirely. 

Because the forbidden region is absent the wave- 
functions will be purely oscillatory. This is similar to 
those wavefunctions obtained when the Wheeler-DeWitt 
(WDW) equation corresponding to a classical signature 
change was solved [Q. It was further found that from 
notions of regularity only a boundary condition analo- 
gous to the Tunneling boundary condition was allowed. 

This article is organized as follows: in section II, we 
describe the minisuperspace model for the case where 
the spacelike sections have zero, positive or negative con- 
stant curvature and we investigate the general behaviour 
of the wavefunctions in such Universes. In section III, 
we discuss the choice of the initial state by making use of 
regularity requirements . In section IV, we make physical 
predictions using the obtained wavefunctions. In partic- 
ular, we compute the probability of creating a Universe 
with a particular spatial curvature and show that it is 
most likely to be an open Universe. The probability for 
inflation to occur is further obtained for any curvature. 
We finish with the conclusions presented in section V. 



II. DESCRIPTION OF THE MODEL 

We consider the quantization of the following FLRW 
metric: 



ds 2 = -N 2 (t)dt 2 



-a 2 (t) 



dr 2 



1 - kr 2 



+ r 2 (d6 2 + sm 2 6d(p 2 )). (1) 



It is well-known that the metric does not fix the global 
topology [jllj. In this article, we will generally assume 
that the spacelike sections are compact (i.e. they have 
a finite volume) with constant curvature characterized 
by k = 0, ±1. Different topologies correspond to differ- 
ent ranges of variation for the coordinates (r, 9, ip). The 
volume of the spacelike hypersurfaces is formally given 
by: 



Vk 



d 3 xVh, 



(2) 



where we have written the metric as: ds 2 — —N 2 (t)dt 2 + 
a 2 {t)h i jdx' l dxi . Then, the Einstein-Hilbert action plus 
the boundary term for this minisuperspace is given by 
the following expression: 



e-h 



, , , , 6fe 6 , a , 
^ fc /d^ ---(-) 



(3) 



The matter is described by a scalar field whose action 
can be written as: 



In order to pass to the Hamiltonian formalism, conjugate 
momenta must be calculated. They are given by: 



c 3 Vk 12ad 
16ttG TV ' 



a 3 Vk <j> 
c N' 



(5) 



^From the last expressions, the canonical Hamiltonian 
can be deduced. It reads: 



H r = N 



I6ttG Ti-2 
c 3 Vk 24a 



C 7I> 

v k 2a 3 WttG' 



We are now in a position where the quantization a la 
Dirac can be carried out cf . [0Jf^] . It consists in replacing 
the two momenta according to the rule: 



d 2 



(7) 



where p takes into account the factor ordering ambiguity. 
The second rule is that the action of the operator H c on 
the wave function ty(a,4>) gives zero. This leads to the 
WDW equation: 

d 2 p d 6 d 2 

q^2 # + 7— *(«. 4>) - — — *(«> 



a da 



Ka 2 dtp 2 



3 ^ ; a 2 (^) 2 [k-(^) 2 ^(a,cp)=0, (8) 



K 2 h 2 C 2 



a 



a Q 



where k = 8ttG/c 4 and do = [kV"(<^)/3] x I 2 . This equa- 
tion is not exactly soluble in its present form so we will 
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make a further assumption. Although we still let the 
wave function keep a scalar field dependence we will ig- 
nore the 2nd derivative term w.r.t <fi in the WDW equa- 
tion. This should be valid during any "slow roll" regime 
where the scalar potential plays the role of an effective 
cosmological constant, i.e. ao = (3/A) 1 / 2 . As a con- 
sequence we find that al/(K 2 h 2 c 2 ) = 3ppi/(5127t 3 ^ p1 pa) 
where lp\ is the Planck length and ppi the Planck energy 
density, ppi = c 7 /(hG 2 ). pa is the energy density of the 
effective cosmological constant, p\ = A/k. Finally, it 
is convenient to work with a dimensionless scale factor 
expressed in units of the Planck length. Therefore we re- 
define the scale factor according to: a — ► lp\a. Equation 
(0) can then be rewritten as: 



d 2 ^(a) 
da 2 



p d*(a) 
a da 

27w 2 pp\ a 2 



-(-)^[fc-(-) 2 Ma) = 0. 
1287T 5 pa o-q ao 

The last equation determines the superpotential: 



U(a;k)- 27U ' 



Ppi / a 
128?r 3 p A ao 



) 2 [fc-(-) 2 



a 



(9) 



(10) 



A related WDW equation has been obtained by using a 
normalization of the scale factor to remove any volume 
factor divergences p2[ . In the following figure, the su- 
perpotentials for the three cases k = 0, ±1 are displayed. 
Note that for k = 0,-1, the superpotentials are always 
negative. This is a crucial difference in comparison with 
the k = 1 case . This means that when k = 0,-1 there 
is no possibility of tunneling anymore since a zero energy 
system is always above the superpotential. This will have 
important consequences which are now investigated. 



solutions. Introducing the variable z(a; k) defined by cf. 



5 (a;fc)-(f) 4/3 ^ /3 (^) 2/3 [fe-(f) 2 ], (11) 

87T pa ao 



the general solution can be expressed in terms of Airy 
functions of first and second kind [fill : 



*(a;fc) 



aAi[z(a; k)} + /3Bi[z(a; k)] _ N(a;k) 
aAi[z(0;k)}+ PBi[z(0;k)} = D(k) ' 



(12) 



where the coefficients a and (3 are arbitrary complex 
numbers determined by the choice of a state for the wave 
function of the Universe. 

We can make some comments at this point. The pres- 
ence of the denominator in the previous equation comes 
from the requirement that the wave function be regular 
everywhere in the minisuperspace. Recall that the wave 
function depends on the scalar field and the approxima- 
tion made previously was only adopted for computational 
convenience. Let us therefore re-consider the full WDW 
equation given by formula (|§|). It is convenient to define 
the quantity a = ln(a/ao). Then the WDW equation 
becomes: 



da 2 
27 v^l ppi 
128tt 3 p A 



(p-1 



i 4a [k 



da 

3 2a ]*(a,0) 



6 <9 2 ^(a, 
k d<j> 2 

= 0. 



(13) 



This equation is separable. If we define \&(a, (j>) 
e~ Q f[a)g{4>) and if A is the separation constant then 
the solution can be written as: 



*(a, 



c(X)e^ a f x (a)e' 



§ \<t> 



AX, 



(14) 
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FIG. 1. Superpotentials for different values of k. The full 
line represents the case k = 1, the dotted line the case k = 
and the dashed line the case k = — 1. 

In order to obtain exact solutions of the WDW equa- 
tion we choose to work with the factor ordering given 
by p = — 1. This permits to work with analytical exact 



where c(A) are a priori arbitrary coefficients and the func- 
tion f\{a) satisfies the equation: 



d 2 A(q) 

da 2 



+ A 



(1-p) 2 



128tt 3 



PA 



fx(a) = 0. (15) 



Note that when considering quantum wormholc solutions 
with massless scalar fields the oscillatory divergence can 
be regulated by a proper choice of the coefficients c(A) 
|l7| |l9| . But we will rather consider single wavefunc- 
tion like solutions for which such a scheme is absent [P0[ . 
When the scale factor becomes small, a <C ao, Eq. ( |15| ) 
can be readily solved. The solution is fx (a) w e lu >-.v a 
where w^,p can be expressed as: 



A 2 



(1-p) 2 



(16) 



In particular when p = — 1, i.e. the factor ordering 
adopted in this article, one has u)\, p =-\ = VA 2 — 1. We 
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are now in a position where the regularity of the wave- 
function in the vicinity of a = can be studied. It cru- 
cially depends on the value of the separation constant, 
see also Ref. [M where a similar treatment has been per- 
formed. If A g] — 1, 1[ then the wavefunction is regular as 
the scale factor goes to zero whereas otherwise the wave- 
function exhibits rapid oscillations. For these values of A 
the wavefunction remains finite. Nevertheless, the "wig- 
gliness" of the wavefunction is also a kind of singularity 
and it seems reasonable not to consider this possibility . 
In this paper, we restrict ourselves to the case A = for 
which the wavefunction is regular at a = 0. This means 
that if? (a — 0, 4>) does not depend on tj> and it justifies the 
presence of the denominator D(k) in Eq. Jl^), see also 
Refs. jjj. It is clear from the previous discussion that this 
is not the most general case since there exists non van- 
ishing values of A such that the wavefunction is regular 
at the origin. However, to our knowledge, there does not 
exist a general study of the boundary conditions for an 
arbitrary value of the separation constant. Moreover, it 
turns out that the two most widely discussed choices of 
boundary conditions, i.e. Hartle-Hawking and Vilenkin 
states, picked out the value A = as will be demonstrated 
below, see also Ref. p0] (in these two cases, D(k) is re- 
sponsible for the appearance of the factors e ±1//v ^ in 
the probability density functions). So as a first approach 
it seems reasonable to consider the case A = only. Let 
us emphasize that Eq. ( [l2|) allows us to study all the 
boundary conditions such that if! (a = 0,4>) is indepen- 
dent of the scalar field and not only the Hartle-Hawking 
and Vilenkin states. 

We should distinguish between the sort of singular- 
ity discussed above, present for example with a massless 
scalar field case , which would tend to be displayed by 
a rapid oscillation in the wavefunction as a — ► (recall ki- 
netic energy ~ "wiggliness" of wavefunction, see eg. 21 1) 
if not regulated and those considered by Ref. [^2] which 
are related to the choice of the factor ordering. There, for 
the closed case, it was claimed that only the no-boundary 
state is regular in the limit a — > for any choice of p. It 
was also shown that the Tunneling solution is regular in 
this limit only if p < 1. This behaviour was caused by 
the wavefunction having growing and decaying, actually 
Modified Bessel functions, solutions in the forbidden re- 
gion. Typically the required Tunneling solution diverges 
like if? ~ a( 1-p '/ 2 in the limit a — > 0. However, this sort 
of divergence is present universally. Solving Eq. (||) for 
k = A = gives the solution 



if! {a) = a +c 2 a J 



(17) 



with c\ and ci arbitrary constants. So the factor order- 
ing divergence remains even in flat empty space. These 
divergences appear intrinsic to these models and one can 
conceive that they should be removed by a renormaliza- 
tion scheme. It was further suggested, in the context of 
quantum wormhole that such factor ordering divergences 
are of no great concern p3[. In practice, as a — > 0, one 



should anyway include a more realistic strong-energy sat- 
isfying matter component. As emphasized by Gott and Li 
Q such a matter source should be expected if only of a 
size due to quantum " zero point" fluctuations. One can 
think of this as being because the scalar potential has 
a "fuzziness" due to quantum uncertainty in the limit 
a — > 0. For the case of radiation, given here by a pa- 
rameter A the WDW equation is then given by, see e.g. 
P5| 



d 2 *(a) pd$(a) 



da 2 a da 
and now with solution. 



A*(a) = 



(18) 



if? (a) = cia^Jp-i (VAa) + c 2 a~^Y p -i (y/Aa), (19) 



where J( p -ru2 an d Y(p-i)/2 are Bessel functions of order 
(p - l)/2 §J. But as pointed out by f|| the Y Bessel 
function still diverges for p > 1 and this term is part of 
the Tunneling boundary condition's solution. However, 
the presence now of oscillatory behaviour complicates the 
adoption of the HH boundary condition. It sometimes is 
given by the prescription "outgoing" plus "ingoing" , if 
so it would also suffer this same divergence, although 
if it was simply the J Bessel function term it would be 
normalizable p6| . But regardless of such considerations 
for the chosen p — —1 case the no-boundary state as well 
as the Tunneling state are immune from factor ordering 
divergences and are a priori both allowed. 

Let us now turn to the study of the behaviour of the 
wave functions given in Eq. (|l2|). The case k = 1 is 
standard. We choose the spacelike sections to be spheres 
although of course other possibilities are allowed. This 
means that v\ = 2-7T 2 . The evolution of the Universe can 
be viewed as the motion of a fictitious particle with zero 
energy in the potential given by Eq. (|1C)|). The particle 
starts to the left of the potential and can proceed to tun- 
nel through the barrier. While in the region a < ao the 
wavefunction exponentially decays. Then, when the scale 
factor is beyond the barrier a > ao, the wavefunction be- 
comes oscillatory. This behaviour is illustrated in the 
following figuren for the boundary conditions a = 1 and 
(3 — i, which correspond to the Tunneling wavefunction 
. The value of p\ is chosen such that ppi/pa = 1000. 
Since we have a 2 , = (3/87r)ppi/ p\, this corresponds to 
a dimensional ao equal to ps 10.9Zpi. This value can be 
thought of as being the size at which the Universe is first 
created. 



*In this article, the figures have been obtained using the 
Mathematica (version 3.0) and IDL software packages 
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FIG. 2. Real part of the Vilenkin wavefunction for k — 1 
and vx = 27r 2 . The value ppi/pA = 10 has been chosen rather 
than the more realistic value ppi/pA = 1000 only for the sake 
of illustration. 

The cases k = 0,-1 are very different. The fictitious 
particle is now always above the potential which is now 
negative. Tunneling is no longer required as classical evo- 
lution is possible. As a consequence the wave function 
always exhibits oscillatory behaviour. To go further, we 
need to know the topology of the spacelike sections. The 
smallest three-hyperbolic manifold is not known. The 
two smallest spaces known are the Weeks space and 
the Thurston space |^8| . The volume of the first one is 
« 0.94 and of the second one is ~ 0.98. The results 
presented here do not depend crucially on the volume of 
the spacelike sections provided they are of the same or- 
der of magnitude. Therefore we choose to work with the 
Weeks space and consequently we take w_i « 0.94. In 
the flat case the volume vq is arbitrary. For definiteness, 
we choose vq = 1. The two following figures show the 
wavefunction for the cases k = and k = — 1 respec- 
tively. The boundary conditions chosen are the same as 
previously. It can be seen that the wavefunctions indeed 
oscillate even before the scale factor reaches the value ag. 
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FIG. 3. Real part of the Vilenkin wavefunction for k — 0, 
vq — 1 and ppi/pA = 1000. 



It can also be noticed that the oscillations in the region 
a < clq are more pronounced for the k = — 1 case than 
for when k = 0. This is due to the different shapes of 
the superpotential. The particle can roll faster down the 
steeper k = — 1 superpotential compared to the k = 
case. 









A A 








V 


J » 





0.0 0.5 1.0 1.5 2.0 

a/a o 



FIG. 4. Real part of the Vilenkin wavefunction for k = — 1, 
v-i » 0.94 and p P i/p A = 1000. 

The fact that the wavefunctions behave very differently 
for different values of k also has an impact on the choice of 
the quantum state. Studying this question is the purpose 
of the next section. 



III. MEASURE AND THE INITIAL STATE 

To be able to make predictions and to calculate prob- 
abilities, we need a suitable measure. If one chooses a 
surface in the minisuperspace perpendicular to the "a" 
direction, then the component of the current associated 
with the WDW equation through this surface is given by 

©: 

j= t -a?(**d a *-*d a **), (20) 

where, in our case, p = — 1. As is well-known this cur- 
rent is not positive definite since the signature of the 
(mini)-superspace is Lorentzian. However, in the WKB 
approximation, where the wave function can be written 
as \ff ~ Ce , this current becomes positive definite and 
permits the calculation of conditional probabilities p9[ . 

Let us first calculate the behaviour of the numerator 
of Eq. (jlj) when the scale factor becomes large. We 
have lim a ^ +oc , z(a; k) = — oo for any value of k. Using 
the asymptotic behaviour of the Airy functions, we find: 

l im N(a; k) = hf^Lgtl ( ( « + ^[-W+if 

a— >+0O 2y/TT \ i 

_(^_^|HM)]^fV (21) 
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Let us now turn to consider the form of the denominator. 
This time the three cases must be treated separately. Let 
us start with the usual case, i.e. k = +1. The value of 
the function z(a;+l) when the scale factor vanishes is 
given by: 



^(0;+i) = (^-) 4/3 («i) 2/3 (^) 2/3 . 

87r p A 



(22) 



Semi-classical considerations are supposed to be valid 
only if we are in a region where ppi/pA 1. We can 
therefore work with the approximation that z(0; +1) S> 
1. In that case, we obtain: 



D(k = +1) 



JL 

2^ 



XOj+lJJ-^eiWOj+D]^. (23) 



The case k = is rather simple since we have z(Q; 0) = 0. 
Therefore, the denominator can be written as: 



D{k = 0) = Ai(Q)(a + V30), 



(24) 



where Ai(0) = 3" 2 / 3 / r ( 2 /3) ~ 0.35502. Finally, we turn 
to the case k = — 1. Now we have: 



z(0;-l) = -(f ) 4 / 3 ( v _0 2/3 (^) 2/3 . 

87T PA 



(25) 



Therefore, in the semi-classical regime, we have 
|z(0; —1)| ^> 1 and z(0; —1) < 0. The presence of the mi- 
nus sign turns out to be crucial. This time the asymptotic 
expansion of the Airy functions for large negative z must 
be used contrary to the case k = 1 where the asymptotic 
expansion for large positive z has been utilized. Thus, in 
this limit, the product D(k — —l)D*(k = —1) is given 
by: 



D(k = -l)D*(k = -1) 



a- 



2^(l+i 2 V-z(0;-l) 
+ l + /7 2 + 2 j ocos-0 ], 



(1 + p 2 - 2pcosip)t 2 + 2(1 - p 2 )t 



(26) 



where t = tan{(2/3)[-z(0; -1)] 3/2 } and (3/a = pe^. 
There is now a danger of obtaining a divergence in the 
current (see the following section) when the product 
D(k — —l)D*(k = — 1) vanishes. The discriminant of 
the second order polynomial in t in Eq. (|2(|) is given by: 



A = -16p 2 sin 2 ip. 



(27) 



Therefore, the polynomial has no real roots except when 
A = 0. This corresponds to the following cases: 



p = 0, and/or ip 



(28) 



where n is an integer. Each wavefunction in the min- 
isuperspace is characterized by the numbers (p, ip) and 
thus belongs to a two-dimensional space. It has been 
shown by Gibbons and Grishchuk |3(]] that this space of 



the wavef unctions is in fact a sphere of unit radius. The 
polar coordinates (0, ip) of a wave function are calculated 
according to the formulas: 



9 = 2 tan" 



ip = ip. 



(29) 



Therefore, the subspace of singular wavefunctions defined 
by Eqns. ( |28| ) is just the great circle going through the 
north and south poles such that ip = 0, n. This shows 
that almost all the wavefunctions are regular except those 
belonging to this subspace. In the case k = 1, the situ- 
ation is different since all the wavefunctions are regular. 
The Tunneling (Vilcnkin's) wave function is such that: 



a = 1, 



/ 8 = i, 



(30) 



which means that its coordinates are (7r/2,7r/2): it is 
regular. On the other hand, the Hartle-Hawking state is 
given by: 



a = 1, 



(3 = 0, 



(31) 



i.e. it is represented by the south pole of the sphere and 
thus belongs to the subspace of the singular wavefunc- 
tions. The divergence occurs when t — — 1 (there is only 
one solution since A = 0) , that is to say when p\ satisfies 
the following expression: 



pp\ 
PA 



8tt 3 

V-l 



32tt j 
3w_i 



(32) 



where m is an integer. This problem is more serious that 
it might first appear because one should not think of 
p\ as being strictly constant (as we have approximated 
during the previous calculations) . It rather slowly varies 
in time during the "slow roll" inflationary period. This 
means that p\ could easily pass through one of these 
dangerous values and so cause a divergence in the wave 
function. Note that this would be across any chosen a = 
const surface, where the semi-classical analysis should be 
valid, and so more serious that the previously mentioned 
factor-ordering type divergences that can occur as a — > 0. 
It is also in addition to the problems of simply having 
real wavefunctions, that do not allow the current j to 
be interpreted directly E^]. This divergent behaviour 
is like that found in the context of the quantization of 
spacetimes which admit a classical change of signature 

To know if a divergent value is actually passed through 
requires a knowledge of the volume For m = 0, 

we have ppi/pa ~ 248/w_i. Depending on the value of 
V—i it could turn out that we are no longer in a regime 
where the semi-classical approximation is valid and so 
the danger could be ignored as being outside the range 
of validity. But for bigger values of m, we certainly will be 
within this region of semi-classical validity. For example, 
if V-i « 0.94 and m = 3, one has ppi/pA ~ 1316. 
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Therefore, it seems that the Tunneling wave function 
can easily be generalized to describe the quantum cre- 
ation of hyperbolic Universes with compact spacelike sec- 
tion whereas the Hartle Hawking wavefunction leads to 
important difficulties, ft is interesting to see that con- 
sidering the quantum creation of compact Universes can 
lead to some progress in the debate Tunneling vs no 
boundary. In the following section, we will consider that 
the wavefunction is placed in the Tunneling or Vilenkin's 
quantum state. 

But before ending this section, we would like to discuss 
the important case of the Hartle-Hawking wavefunction 
in more detail. We will concentrate on the case k = — 1 
versus k = 1. The first step of the derivation consists 
in establishing the Euclidean equations of motion. With 
the help of Eqns. (||) and (^), it is easy to see that in the 
gauge N = they read: 



a 
a 



N 2 



kV ((/>), 



+ 3— (j)' - N 



3 



0. 



6 



k-, + —v(4>) 



o. 



(33) 
(34) 

(35) 



where a prime denotes derivative with respect to t = —it. 
The hypersurface on which we evaluate the wavefunction 
is the hypersurface such that r = 1 for which we have: 
a(r = 1) = a and (f>(r = 1) = <fi. The no boundary 
conditions are [15]: 



a(r = 0) = 0, 



d0 

7w 



(r = 0) = 0. 



(36) 



In the slow roll regime, the solutions satisfying all the 
boundary conditions are found by integrating the Eu- 
clidean equations of motion. They are given by: 



a(r) 



sin N 



sin N 



(r)=0. (37) 



The value of k does not appear explicitly in these equa- 
tions but is in fact hidden in the algebraic equation sat- 
isfied by the lapse function: 



N 



na 2 V 
3k 



(38) 



[15]. Everything is now present and the Euclidean action 
can be calculated along this solution. Then the Hartle- 
Hawking wavefunction is given by $hh(o, 0) ~ exp(— 5e) 
(we have dropped the bars in order to avoid cumbersome 
notation) and can be expressed as: 



*HH(a, 4>) ~ exp 



6i>i 
n 2 cV 



, | | 2 KV 

1 — 1 — a — 



3/2- 



(40) 



In the region where a <^ ao the previous wavefunction 
can be written as \l/HH(a, (j>) ~ exp[(3?;ia 2 )/(K;c)] and we 
see that this is indeed independent of the scalar field and 
therefore that we have picked out the value A = 0. 

Let us now turn to the case k = — 1. The solution to 
the algebraic equation giving the lapse function can be 
written as: 



sinh 



(41) 



Interestingly enough the solution is now unique. The ex- 
plicit expression of the scale factor can be easily deduced. 
It reads: 



(42) 



«( T ) = V ^7 sinh 



sinh a\ 



In the same manner, we can evaluate the Euclidean ac- 
tion and therefore find the corresponding no boundary 
state. One obtains: 



*HH(a,<?!>) 



6v_i 
6XP i l ^cV 



,kV 



3/2-1 



(43) 



where we have again suppressed the bars. In the limit 
where the scale factor goes to zero, this wavefunction 
can be written as 'I'hh^, 0) ~ exp[— (3iv-ia 2 ) / (kc)]. It 
confirms that we have picked out a vanishing separation 
constant. This also shows that in the Euclidean region, 
the wavefunction is oscillatory. One could easily find the 
no boundary wavefunction in the classical region with the 
WKB matching method. Of course, we would recover the 
fact that the denominator is now a trigonometric function 
which was the main reason for arguing that the Vilenkin 
state is preferable. In fact we can go a step further and 
say that the very concept of a no boundary wavefunction 
is in danger by this kind of analysis. Indeed, from Eq. 
( |35| ) we see that the no boundary conditions (|36|) imply 
that: 



At this point, one should consider the two cases sepa- 
rately. The case k = 1 is standard. The solution of the 
previous equation can be expressed as: 



1. _ x ( IkV 

(n + — J7T ± COS I (/ 



(39) 



It is common to consider only the case n = together 
with the minus sign, for a fuller discussion see Ref. 



N 2 { dr J 



(44) 



N 2 dr 2 



2 (r)d«l to be 



In order for the metric ds 2 
Euclidean and regular at r = we need a(r) sa Nt 
which implies k = +1! In the case fc = — 1, it is easy to 
see that we have in fact a(r) sa —iNt. This means that 



ds^ 



N 2 (dT 2 



2 dil 2 ). The corresponding manifold is 



no longer Euclidean. In addition, the metric found above 



7 



is very reminiscent of that used by Hawking and Turok 
[9] in their analysis. Therefore it seems that there exists 
an interesting link between the approach advocated here 
and the one of Ref. [9]. This requires further studies 
which are beyond the scope of the present article. 

To end this section, let us come back to the question 
of the regularity of the wave function. We have just 
shown that the requirement |\I/(a;fc)| < oo everywhere 
in the minisuperspace and for any k favours the Tun- 
neling wavefunction over the Hartle Hawking state. On 
the other hand requiring that the wavefunction be regu- 
lar as a — ► for every factor ordering leads to a Hartle 
Hawking state for the k = 1 case p2[ . Therefore, the two 
requirements are strictly not compatible if the analysis 
was extended to p > 1. However, we have previously 
mentioned that in more realistic models the actual pres- 
ence of non-inflationary matter could also complicate the 
adoption of the HH wavefunction. 

Regardless of these complications, and because factor 
ordering problems, which occur in the unrealistic limit 
a — > 0, appear less serious than the divergences found in 
this section for k = — 1 models, we will next consider the 
predictions with these Tunneling wavefunctions. 



IV. PREDICTIONS 

In this section, we address the problem of computing 
physical predictions from the previously obtained wave- 
functions. It is well known that, due to the fact of ob- 
taining non-normalizable wavefunctions, this is a difficult 
task. In Ref. , a method to overcome this problem has 
been proposed. The idea is to use the current defined by 
Eq. (gg) in the WKB regime. This leads to well-defined 
probabilities. Using the asymptotic form of the wave- 
function given by Eq. (|2l]), the following expression for 
the current, valid for any k, can be found: 



2 3. 
7r 8n' 



1/3. 



2/3 



D(k)D*(k) y PA 



(PPl)-l/3 



(45) 



As expected this expression no longer depends on the 
scale factor. Going further requires a knowledge of D(k) 
so we must now treat the three cases separately. Let us 
first consider the case k = 1. The expression of D(k = 1) 
can be deduced from Eq. This leads to: 



-1 



3i>i 3 "i ppi 

e 16ir 2 PA . 

7T 



(46) 



We recover the well-known expression for the closed case. 
The coefficient of proportionality is not of great interest 
since the probability distribution is not normalizable and 
instead must be used to calculate conditional probabil- 
ities, see below. The ratio ppi//9A only appears in the 
argument of the exponential function. This is because 
the factor [z(0; + 1)]~ 1/2 in the D(k = l)D*(k = 1) term 



cancels exactly the term (ppi/pj 



v-l/3 



in Eq. fl45h. The 



expression of the current in the case k = can be estab- 
lished from Eq. (|24|). It reads: 

^=»>=w<!> 1,a »» ,3 ,/3 - (47) 

Since D(k = 0) is just a constant, this does change the 
dependence in (pp\/ Pa)' 1 ^ 3 of Eq. (M). Finally, we 
turn to the case k = — 1. From Eq. ((26|), we have that 
D(k = -l)D*(k = -1) = (1/tt)[-z(0;-1)}-^ 2 . This 
time no exponential function appears as was the case for 
k = 1. This is due to the fact that we have chosen the 
boundary condition that D(k = — 1) simply be a phase. 
However, now for the case k = — 1, the factor Ppi/pa 
simply cancels out. As a consequence, one finds: 



-1) 



3«_i 

47T 



(48) 



^sup 

is the value 



The current turns out to be independent of p\. 

As already mentioned, the current can now be used 
to calculate conditional probabilities. For example, the 
probability of having an initial value <f>{ of the scalar field, 
greater than the value needed to solve the problems of 
standard cosmology, <^ su f, knowing that < (f> < 
can be allowed. In the last inequality 4> sup 
at which semi-classical considerations cease to be valid, 
i.e. when the potential reaches the Planck scale, ppi = 
p\ = V(4>) = m P1 in the Planck system of units. Also, 
the fact that the minimal value for <fi is zero is not a prob- 
lem here because the wave function is the Tunneling one. 
This would no longer be true if the state were the HH 
one. We will evaluate the conditional probabilities for 
the prototype chaotic inflationary scenario, with scalar 
potential of the form V(<j>) = (A/4!)</> 4 . In this context, 
the initial value of t he field necessary to get N e- folds 
is (f>i = \J (N + l)/7rmpi. If we consider that sufficient 
inflation is obtained when N = 60 then we find that the 
scalar field has to start at <t>i = cW = 4.4mpi. The value 



of <^ sup is given by 



'sup 



: (fsuf 

24V4 A -i/4 mpb For c haotic in- 
flation, one has A ~ 10~ 15 in order to reproduce the value 
of the quadrupolc of the Cosmic Microwave Background 
Radiation anisotropy, Qrms-ps, measured by the COBE 
satellite, i.e. Q rros _ps w 18 x 1CT 6 K This implies 
0sup ~ 1.2 x 10 4 topi. The probability to have sufficient 
inflation in these three cases is computed according to 
the formula |15|: 



P{k; <j>i > <t> But \0 < <k < <t> Bup ) = P{k) 



(49) 



It will be more convenient in the following to rewrite 
the previous definition of the conditional probabilities as 
P(k) = 1 - R(k; N, A) where R(k; N, A) is given by: 



R(k; N, A) 



(50) 
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We are now going to calculate R in the three cases. 
Let us start with the case k = 1. The current 
is given by Eq. (fl6|). Using the change of variable 
u = %rr4J(\cj) 4 ) and the formula (3.381.6) of Ref. || 
one can easily show that the denominator is given by 
(24/9) 5 / 8 e- 3 / 16 W_5 _i(9/24) where W^ v is a Whit- 
taker function. Since W_&_ 1 (9/24) pa 0.52 we find that 
the denominator is approximately equal to 0.80. In the 
same manner, the numerator is equal to: 



9tt 2 



-i -5/8 



A(^V + 1) 2 



( 9ir 2 



(51) 



Since the parameter A appears at the denominator of the 
argument of the Whittaker function, this has a very large 
value. Therefore, we can use the asymptotic expansion 
of the Whittaker function for large values of its argument 
given by |1: 



lim W lt , v (z) = e-iz^[l + 0{-) 

\z\— »oo \ Z 



(52) 



Thus, the value of the function R(l; N, A) can be ex- 
pressed as: 



R(1;N,X) pa 0.0045 \*(N + l) 5 e i <«+ 1 ) 2 



(53) 



This expression is valid for any value of N and A pro- 
vided that A is a small number. Putting N = 60 and 
A = 10~ 15 , wc find that R(l; 60, 10~ 15 ) pa 10 



-1Q 1 



an 



extremely small number which has its origin in the pres- 
ence of the parameter A in the exponential factor. We 
conclude that P(l) pa 1. The calculation of P(0) is eas- 
ier. Using Eq. ft?]), we find that: 



R{0;N,X) = 24~* ( ) A^. (54) 



This gives i?(0; 60, 10~ 15 ) pa 1.5 x 10~ 9 leading to P(0) pa 
1 again. Finally, the case k = — 1 is straightforward. 
Using Eq. (ftsj), we can easily establish that: 



R(-1;N,X) = 24^3 



ifN + 1 



A3. 



(55) 



This results in R(-l; 60, 10~ 15 ) pa 1.67 x 10~ 4 . There- 
fore, one can say that P(— 1) pa 1. The conclusion is that, 
in the three cases, the probability turns out to be close to 
one. This means that sufficient inflation is a prediction 
of the Tunneling wavefunction whatever the value of k 
is. In this respect, the three cases are equally compatible 
with there being a near definite prediction of inflation oc- 
curring. This is actually an significant improvement over 
classical notions of whether inflation will occur. Which, 
due to an infinite divergence over an arbitrary scale fac- 
tor, gives an ambiguous prediction even for apparently 
inflationary potentials 0. We note however than the 



probability is closer to one in the case k = 1 than in the 
cases k = 0, —1. 

We can pursue this reasoning a step further by also al- 
lowing k to be quantum (q) variable. Such a possibility 
could be a first step in modeling quantum fluctuations 
in the geometry. Unlike the classical case where the cur- 
vature is simply a fixed (c) number constant we are now 
assuming that the curvature also is in a quantum en- 
semble. The conditional probability for having a given k 
knowing that k = 0, ±1 can be calculated for a fixed value 
of the scale factor. This probability is formally defined 
according to the equation: 



P(a;k) 



l*(a;*)l s 



E I= o,±il*(«;0l s 



(56) 



The use of this definition requires some comments. Let 
us first recall how the definition of Eq. (49) can be jus- 
tified. Let M. be the minisuperspace and .Mwkb £ M. 
the region of the minisuperspace where the wavefunction 
can be well approximated by the WKB wavefunction. In 
non relativistic quantum mechanics, the time component 
of the current, |\&| 2 , gives the probability density func- 
tion. Since the wavefunction is normalizable, the proba- 
bility of finding a particle in the interval [a, b] is given by 
/ dxj*! 2 / da;|\I>| 2 . In quantum cosmology, the cur- 
rent is positive definite only in the WKB regime. In this 
regime, \P(a; k) is not normalizable, i.e. f M dq a \^f\ 2 , 
where dq a is the volume in Ai, is not finite. This does 
not mean that J M n(q a )dq a \^\ 2 is infinite since the mea- 
sure n(q a ) is a priori not known. In the WKB regime, the 
wavefunction is by definition peaked over AIwkb- There- 
fore, one has J M ^q a )dq a \^\ 2 Pa dg Q |*| 2 . Thus, 

the same rule as in ordinary quantum mechanics says 
that the probability of finding the system in the region 
UeMwKB is given by f K dq a \V\ 2 / dq a \^\ 2 . Eq. 

( [ilil ) is a special case of this more general formula. 

The interpretation of Eq. (|5^) is roughly the same. 
Although it is, of course, more contentious to apply this 
reasoning to obtain the geometry, here represented by 
k, compared to how we previously obtained the initial 
matter component <p. The matter calculation is a fairly 
straightforward adaptation, as stated above, of usual 
quantum mechanics reasoning whereas the quantization 
of the geometry might require more extensive alterations 
to quantum mechanics. We will proceed with the no- 
tion that k is not initially fixed but is rather undefined 
in a quantum state with "equipartition" among all pos- 
sible k. In the context of the histories approach of quan- 
tum mechanics, P(a; k) represents the probability that 
the Universe "choose" one of these possible histories. It 
will not come as a surprise that the Universe goes down 
the lowest potential case. 

We would like to emphasize how P(a; k) differs from 
the notion of the probability of a change of topology once 
the Universe has been created, i.e. the probability of hav- 
ing, for example, k = 1 for some value of the scale factor 
and then k = for another (bigger) value. If topology 
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changes are required from one classical model to another 
then an explicit time dependent curvature should be in- 
troduced. Such a construction with dk/dt ^ has been 
explicitly made in Ref. [^3|] and it has been shown that, 
in this case, a passage to a midi-superspace description 
is mandatory. However in the present example we are 
assuming that the curvature is first a quantum variable 
and so not fixed in a particular classical state. What clas- 
sical curvature state the universe first evolves to is our 
present concern, not whether topology changes still occur 
once this classical state is first achieved. As the universe 
becomes increasingly classical, by for example gaining en- 
ergy from falling down the potential, then the curvature 
will no longer be in a quantum superposition but will be- 
come increasingly in a specific classical state. How this 
"measurement" takes place will be a rather complex pro- 
cess and likely to depend on quantum mechanical inter- 
pretational questions. But for our present purposes, this 
rough notion that the curvature will eventually "crystal- 
lize out" into a classical state should suffice. 

On the following figures, the evolution of the three 
probabilities as a function of the scale factor are dis- 
played. The first figure is for the case ppi/pa = 10. 
This means that the dimensional quantity ao is equal to 
ao w lp\. Let us also recall that we have taken v\ — 2ir 2 , 
v o = 1 and V-i « 0.94. The choice of the volume can in- 
fluence the behaviour of the probabilities for small scale 
factors but in the limit of big scale factors, they are 
mainly determined by the ratio ppi//9A- 
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FIG. 6. Probabilities P(k) for p P1 /p A = 100. The solid line 
represents the case k = 1, the dotted line is the k = case 
and dashed line is the k = —1 case. 

Finally, the third and last figure represents the case 
Ppi/pa = 1000 which corresponds to ao ~ 10.9Zpi as al- 
ready mentioned. This case could be considered as the 
most realistic one since inflation takes place at an energy 
comparable to the GUT scale, i.e. pa ~ 10 16 GeV. 
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FIG. 7. Probabilities P(k) for p P1 /p A = 1000. The solid 
line represents the case k = 1, the dotted line is the k — 
case and dashed line is the fc = — 1 case. 



FIG. 5. Probabilities P{k) for p P i/p A = 10. The solid line 
represents the case k = 1, the dotted line is the k = case 
and dashed line is the k — — 1 case. 

The second figure represents the case where ppi/pA = 
100. This corresponds to ao ~ 3.5Zpi. 



Let us now comment on these figures in more detail. 
In each case, we can point out the following features. At 
vanishing scale factor, the probabilities are all assumed 
equal to 1/3 for any value of k. Provided the probabili- 
ties are roughly equal the actual values will not affect the 
predictions significantly. In the region where the scale 
factor is still of order ao, the behaviour of the probabili- 
ties are rapidly evolving with a. Finally, when the scale 
factor becomes large in comparison with ao, the proba- 
bilities tend to a constant value which depends on the 
ratio pp\/ pa- This suggests some interesting ideas. For 
example, the fact that the probabilities remain similar in 
the region a < ao could mean that here the topology can 
easily fluctuate due to quantum effects. On the other 
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hand, when the scale factor is such that a > oj, the 
probabilities differ significantly and are almost constant 
as a function of a. This can be easily understood if one 
looks at the asymptotic behaviour of |\I>(a;A;)| 2 when a 
becomes large. Using Eq. (|l2|), one finds that: 



\V(a;k 
\V(a;k 
|*(a;fc 



4 , 3 



a 87r 
1 



PA 



^1/6 

47rAi 2 (0) V 3 ; 



1(^)1/2(^1/2 

a v 8vr 7 >a ' 



( m }1/ 6 5 



!/3 • /). 



(57) 
(58) 

(59) 



In each case |^(a;/c)| 2 behaves as 1/a and therefore 
P(a; fc) becomes independent of the scale factor when 
a a . Fluctuations in topology are more likely, re- 
gardless of the scale factor, provided the energy density 
is near Planck values. But as the value of p\ reduces this 
rapidly becomes less likely. When pp\/ pa ^> 1, it is clear 
from the figures that a definite prediction can be made 
since one of the probabilities becomes equal to one: 



P(-l)«l, P(0)«0, P(+l) 



0. 



(60) 



The fact that the greater the ratio pp\/ pa is, the better 
the prediction is [i.e. the closer to one P(— 1) is] means 
that when inflation takes place at energy well below the 
Planck scale (typically the GUT scale seems to be the 
most physical case), topology changes become strictly 
forbidden as the scale factor becomes large. In practice 
as the scalar field rolls down the potential the effective 
cosmological constant pa reduces and topology changes 
become increasingly unlikely. 

Finally, let us comment on the prediction itself. It 
says, since an inflationary phase is expected to give a fi- 
nal a ~ 10 3 %i, that our Universe is most likely to be 
open, or at worse flat, a quite interesting statement in- 
deed. Bearing in mind that we have not tried to allow for 
the more numerous hyperbolic topologies corresponding 
to k — —1, we have probably still underestimated this 
overwhelming probability of obtaining an open universe. 
However, it does not say that ilo is necessarily much less 
than unity and so could actually be distinguished from 
fio = 1 by observations. This is because we have seen 
that even in the open case, sufficient inflation is a predic- 
tion of the Tunneling wavefunction. Thus, inflation will 
probably drive f2i n i < 1 to a value very close to one even 
at the present age of the universe. The actual form of 
the scalar potential is crucial for determining such prop- 
erties, a quantum description alone does not mandate 
such values for f^o- Once the inflationary phase finishes 
the matter behaves effectively like that of dust 0] so now 
obeying the strong-energy condition. Eventually the cur- 
vature term will start again to dominate the dynamics, 
with the open model expanding infinitely into the fu- 
ture. Curvature can play an important role at both the 
beginning and end of the universe's evolution because of 



the drastically different change in the behaviour of mat- 
ter: from being a cosmological constant to becoming like 
dust. Any small residual cosmological constant would 
also come to dominate at later times. Recent develop- 
ment, although still contentious, have suggested such a A 
term is necessary to explain the Supernova data |m| , but 
such a value is still extremely small pa ~ 10~ 120 . Under- 
standing such fine detail while quantum cosmology takes 
a rather "broad brush" approach to calculating the vari- 
ous quantities remains a difficulty. We would just remind 
readers that the wormhole, and related mechanisms, that 



were suggested should predict A = exactly 35 j, could 
instead give other values which are only approximately 



zero -see e.g. 



V. CONCLUSIONS 

In general relativity the global topology has to be im- 
posed as an initial condition. In quantum cosmology 
most studies have concentrated on elliptic space which for 
the FLKW model is simply a closed universe. However, 
hyperbolic and flat spaces are mathematically speaking 
more numerous and many are also compact. These give 
rise to geometries that are locally described by the FLRW 
metric with k = —1,0. For the simple DeSitter model the 
k = —1,0 cases are no longer distinguished by the pres- 
ence of a forbidden or Euclidean region at small scale 
factors. 

Because the Euclidean nature of the model is now ab- 
sent it might seem that the smooth geometric picture 
of the Hartle-Hawking no-boundary proposal is a serious 
loss. But a major weakness of "quantum creation of the 
universe" ideas is that the forbidden region is anyway 
eroded by the presence of strong-energy satisfying mat- 
ter. Such matter will generally prevent the Euclidean na- 
ture of the model for small scale factors. But on general 
grounds such matter should be present, if only because of 
"zero-point" fluctuations. A Lorentzian region will tend 
to form anyway as a — > 0. 

With a scalar field source a possible problem still re- 
mains since a singularity caused by the kinetic energy 
of the scalar field blowing up could occur. The kinetic 
energy of the scalar field behaves as a 'stiff' equation of 
state. This causes an effective ~ —a~ 2 term in the WDW 
potential that would tend to create an infinite 'wiggli- 
ness' in the wavefunction as a — > 0. By ensuring that the 
wavefunction be independent of the matter as a — > this 
possible problem is avoided. But one must bear in mind 
that one is removing the singularity by fiat, one should 
not claim that the process of quantization alone is achiev- 
ing this, as many studies erroneously assume. Likewise 
in the usual closed models the Hartle Hawking boundary 
condition only gives a Euclidean region by imposing that 
the matter fields are not allowed to dominate as a — > 0. 

Once it is accepted that a Lorentzian region is any- 
way present for the smallest scale factors it is no longer 
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a major fundamental difference whether one then has a 
forbidden region away from the origin. In the flat and 
open cases the forbidden region indeed goes away and 
we have argued that these cases are in a sense more 
favourable. Further, in the open case the requirement 
that the wavefunction should remain finite at arbitrary a 
allows boundary conditions that include Vilcnkin's "out- 
going only" but not the no boundary choice. 

For the specific Tunneling boundary condition, com- 
parisons between models that only differ by their value of 
k can be made. If the initial curvature is given by a quan- 
tum ensemble, so allowing any possible k, then one can 
conclude that open universes are more likely. Although 
the models considered have only a cosmological constant 
as their matter source. This seem to contradict the no- 
tion that for DeSitter space "all curvatures are equiva- 
lent" . But because we consider the evolution from a = 0+ 
and regularized the wavefunction around this point we 
have derived "propagators" to go to arbitrarily large scale 
factors. The universe is given a choice from its conception 
which path to follow. In the closed model one can think of 
the universe being "held up" waiting to tunnel through 
the barrier. While the open models gain a "push" in 
falling down a steeper WDW potential compared to the 
flat case. Once the scale factor becomes sufficiently large 
the possibility of topology change depends entirely on the 
energy density of the scalar field driving inflation. As this 
reduces below Planck values the prediction rapidly gives 
that an open k = — 1 universe is favoured. Interestingly 
at small scale factors and/or Planck energy densities the 
universe might "pin ball" between various possibilities 
before settling into one final curvature. The "no hair" 
property of DeSitter space, of having finite causal hori- 
zon, might also allow various regions to develop different 
curvatures. Although the open case would still domi- 
nate pockets of closed curvature could also be created. 
This could depend on quantum interpretations, whether 
"many worlds" or "single-history" quantum theories are 
possible cf. |37| . 

Regardless of the actual value of curvature inflation 
is predicted in all cases. This resolves an ambiguity in 
classical measures of inflation. If quantum cosmology did 
nothing else but gave a definite predictions for inflation to 
occur it would be very significant. The only drawback is 
that the scalar potential must be chosen to be of the cor- 
rect shape, just as the potential of the Hydrogen atom 
has to be provided before one does any quantum me- 
chanics. The potential is further constrained by the need 
to create sufficiently small fluctuations and gravitational 
waves: roughly speaking inflation should occur at GUT 
~ 10 14 GeV energy scales 38 3?|. Ultimately inflation is 
a classical phenomena that can not be driven by quan- 
tum notions alone. We have modeled curvature as being 
described by a quantum variable and subject to notion of 
probability analogous to the usual initial matter distri- 
bution calculation. There is the worry that in trying to 
determine the curvature, which is also a part of the po- 
tential, we are going beyond the scope of what quantum 



mechanics should be used to predict. We have assumed 
that the universe is free to take the path of least resis- 
tance i.e. follow the open case, but maybe there is actu- 
ally no freedom in this choice and it is pre-ordained what 
curvature the universe should take before the universe 
comes into existence. This is a rather deep problem that 
affects quantum cosmology in general, what variables are 
free to roam and which are fixed externally imposed con- 
stants ? Is there a classical scaffolding surrounding the 
initial universe or is every variable initially a quantum 
variable? Quantum cosmology seems ambiguous why cer- 
tain variables (e.g. </>) are given by distribution functions, 
although constrained by the chosen boundary conditions, 
while others (e.g. k) are imposed with no longer appar- 
ent quantum uncertainty. Recall that all classical notions 
seem suspect as the Planck epoch is approached cf. [[Io| . 
But in the meantime it seems that, in admittedly sim- 
plistic models, different curvatures can be considered and 
an argument made that the open case is favoured. If the 
curvature is initially fixed and so not subject to quan- 
tum uncertainty, then one can still argue that open uni- 
verses are as possible as closed ones with Tunneling like 
boundary conditions. Understanding better the measure 
of possible topologies would seem the next helpful step to 
see if arguments can be made to favour a specific choice 
cf. @ 

In summary, we have considered quantum cosmologi- 
cal models with arbitrary curvature. Although the mod- 
els are all compact they are all possible candidates for 
quantum creation of the universe models. Unlike always 
starting from closed models which have a forbidden re- 
gion at small scale factors, one can work directly with 
the curvature of one's choice. Compared to recent in- 
stanton methods of creating an open universe from an 
initial closed one, one cuts out the unnecessary closed 
stage. One has the further prediction that open universes 
are favoured followed by flat ones provided in some sense 
the universe has the choice of deciding its curvature. The 
closed universe case is strongly suppressed in comparison. 
One can still obtain f2o ~ 1 m &U cases since a long pe- 
riod of inflation is strongly predicted. Unfortunately this 
large inflationary period would appear to wipe out any 
interesting "multiple images" due to topological effects 
that are presently being searched for. Only for lesser 
amount of inflation would such effects be apparent in 
the patterns of Cosmic Background radiation or multiple 
galaxy images -see 111]]. But interestingly with an infla- 
tionary matter source, the curvature still dominates the 
dynamics both at the beginning and the end of the uni- 
verse, this seems reasonable that the end of the universe 
should reflect its origins. This is unlike the conventional 
big bang model where the curvature only dominates in 
the far future and "only matter matters" during its initial 
phase. 
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